Entanglement between two fermionic atoms inside a cylindrical harmonic trap 
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We investigate quantum entanglement between two (spin-1/2) fermions inside a cylindrical har- 
monic trap, making use of the von Neumann entropy for the reduced single particle density matrix 
as the pure state entanglement measure. We explore the dependence of pair entanglement on the 
geometry and strength of the trap and on the strength of the pairing interaction over the complete 
range of the effective BCS to BEC crossover. Our result elucidates an interesting connection between 
our model system of two fermions and that of two interacting bosons. 
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I. INTRODUCTION 

Recent experiments with lattice fermions across Fesh- 
bach resonance have raised significant hope for the ap- 
plication of atomic quantum gases to implementations of 
quantum information processing [H, 0, Hj] . These experi- 
ments usually begin with the preparation of two fermions 
into each optical lattice site. Making use of a Feshbach 
resonance [4(, the two-atom scattering length is tuned 
from a small positive value to the attractive side by vary- 
ing an external magnetic field. During this process, the 
motional states of the fermionic atom pair in each lattice 
site can occupy different Bloch bands, exhibiting new fea- 
tures that need further theoretical explanation Q. 

In quantum information science, quantum entangle- 
ment is viewed as the enabling resource for its poten- 
tial power over classical information processing, as in 
quantum teleportation [f| 0, H, H, EH and cryptography 
[ill [H, EH- In this paper, we present a thorough in- 
vestigation of the pair entanglement of two fermions in 
a single optical lattice site, approximated as a cylindri- 
cal harmonic trap. Similar studies have been carried out 
in a one dimensional harmonic trap 1141 and in a three 
dimensional spherical harmonic trap [151 ]. The present 
study differs in two significant aspects: first, we extend 
earlier studies to the case of three dimensional cylindri- 
cal harmonic traps; second, we allow for the possibility 
that the two interacting atoms may form a molecular 
bound state, albeit in the broad resonance regime. In 
the model to be presented below we will parameterize 
the two atom interaction using a general formalism com- 
monly adopted in the many body system of BCS to BEC 
crossover 

(3 m h m. 

For a broad Feshbach reso- 
nance, where most of the currently employed resonance 
theories do fall into, the earlier results based on a sin- 
gle channel model are easily recovered by excluding the 
molecular component. 

This paper is organized as follows. First, we present 
the model system and our formulation, essentially paral- 
lel to the development and notations of Ref. ||. Next, 
we focus on the broad resonance regime and briefly dis- 



cuss the molecular component. This paves the way for 
the discussion of our central result — a thorough inves- 
tigation of the pair entanglement and its dependence on 
trap strength and geometry, and on atom-atom interac- 
tion strength. Finally, we conclude with a brief summary. 



II. TWO INTERACTING FERMIONIC ATOMS 
IN A HARMONIC TRAP 

Following the successful model proposed in Ref. [f|, 
two fermions are assumed to be located in the same lat- 
tice site, which for this study is approximated as a har- 
monic trap. Each optical lattice site represents an inde- 
pendent system, since the quantum tunnelling, or hop- 
ping between neighboring sites, is negligible for deep op- 
tical lattice potential of interest. The Hamiltonian for 
the system of two spin-1/2 fermionic atoms is given by 



H 



Ema^ama- + vb^b 

mcr 



(1) 



where E m = J2j= x y z^ji 171 ] + V^) is the harmonic 



oscillator energy for state labelled by m 



,m z ) 

with angular frequencies (ui x , oj v ,u) z ). a) ma is the creation 
operator for a fermionic atom in the open channel with 
energy E m and spin a. 6 T is the creation operator for the 
two-atom bound state, a bosonic molecule in the closed 
channel with its center of mass wave function fixed ex- 
actly at the harmonic ground state, a result of the simple 
approximation that the optical lattice trap potential for 
the bound state molecule is simply equal to the sum of 
the trap potentials of the two atoms, v is the energy 
difference between the closed channel bosonic molecule 
and the two fermions in the open channel. a m>n is the 
coherent coupling element converting two open channel 
fermions into a closed channel bosonic molecule, which is 
defined to contain a common constant prefactor a, i.e., 
c*m.n = a(O c .m. , V'rcil 111 ! n ) with |m) being the harmonic 
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orbital of a single atom. The relative motional part of 
the molecular ground state, |^ re l)j will be approximated 
as a contact 6(r) function, since it is typically of atomic 
size [jl HFj], much less than other length scales in this 
problem. Such a simplification contains an ultra-violet 
divergence that can be removed by a suitable momentum 
cutoff with a renormalized detuning [20| . An alternative 
formulation involves the use of the regularized delta func- 
tion [2l|. The cylindrical harmonic trap is characterized 
by the trap frequencies u x — uj y = lu±_ = lo z /X with A 
parameterizing the trap aspect ratio. Within this model, 
the two fermionic atoms in the open channel (of being 
atoms) at bands m and n can be converted into a closed 
channel bosonic molecule, or vice versa. The eigenstate 
of Hamiltonian ^ can be expressed as 



I*) = ( P b * + £ ^ a U a ll I i vac > 



(2) 



again following Ref. The various coefficients are de- 
termined from the following coupled equations 



??m,n 

E —V 



= 



a 



m,n 



E 



E-E; 

a 
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E-EL 



1 \ - 

7? = 1 + ^ 



(E — £m,n) 



(3) 
(4) 

(5) 



As in Ref. the divergence in Eq. (TJ| can be re- 

moved by introducing a cutoff in the summation. Such 
a procedure will renormalize v to v* . By comparing the 
free space expression of the energy E (renormalized) to 
the results of resonance scattering at low energies, the 
parameters v* and a are then matched to the experi- 
mentally relevant parameters a s and r , where a s is the 
s-wave scattering length between fermionic atoms in dif- 
ferent internal states and ro is the effective range. The 
result is @ 



v*\ro\m 
2ft 2 ' 



(6) 



where m is the mass of the atom. A dimensionless pa- 
rameter x — E/(2ftuj±) — 1 — A/2 is then introduced, with 
which the energy quantization condition becomes 



d ±- | M 

a s d± 



1 



A 



:F 



X 1 

A' A 



,(7) 



where d± = \fhj (mu}±), and the function on the right 
hand side of Eq. is defined in Ref. [22] 



F(u,tj) 



dt 



■qe 



VI -e-*(l 



-¥ 



1 



(8) 



For a spherical trap this reduces to the well-known result 
of F(-x,l) = -2^T(-x)/r(-x- 1/2) [23]. Further- 
more, if ro = 0, the energy spectrum and eigenfunctions 



coincide with the results of two atoms in a harmonic trap 
as studied previously in Refs. [13,[2^|. This is a straight- 
forward conclusion, since for ro = 0, the current model 
describes the same ph ysic al process quantified by a single 
s-wave parameter a s [22L |23| . A non-zero ro as incorpo- 
rated in Ref. [5] allows for a more general model in- 
cluding both open channel fcrmions and a closed channel 
bosonic molecule. 
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FIG. 1: (Color online) Energy spectrum of two fermions in 
a cylindrical harmonic trap versus inverse scattering length 
in the broad resonance regime. Here we choose A = 5/6 and 
|ro|/d± =0.04. 

The energy spectrum versus inverse scattering length 
— dj_/a s for the current model is shown in Fig. [TJ Across 
a Feshbach resonance, the scattering length can be tuned 
according to a s = a^ & {l — A / (B — B )) [2J|, where a^ g is 
the background scattering length. Bq and A are the res- 
onant field and width, respectively. The effective range 
ro is found to satisfy ro = — 2ft 2 /(m/j,ab g A) [25|]. /i is 
the magnetic moment difference in the open and closed 
channel. 

In this work when dealing with a broad Feshbach res- 
onance, the validity of our model calculation requires 
that I ro I <C d±, \a s \. As in previous studies [22|, HU, the 
eigenenergy remains an increasing function of the inverse 
scattering length —d±/a s . In a more complete treatment, 
the closed molecular channel is observed to differ from 
the lowest two atom (bound) state |e-a) as shown in this 
figure in the internal spin state ■ 



III. RESULTS ON PAIR ENTANGLEMENT 

In this study, we will limit our discussion to the pair en- 
tanglement for a broad Feshbach resonance, correspond- 
ing to the regime of |ro| -C dj_< The nonclassical cor- 
relations (squeezing and entanglement) were previously 
considered for non-condensate atoms across a Feshbach 
resonance in free space [26]. In the following we shall 
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therefore focus on the two adiabatic eigenstates labelled 
as |a-d-b) and | b-f-c) in Fig. [TJ 

First, we briefly review the result on the molecular 
component according to Ref. Starting from a small 
positive scattering length, when adiabatically following 
the state |a-d-b), the molecular component P mo i = 
1 remains small in the broad resonance regime 
More precisely, 



f3 2 



dx 



(9) 



i.e., the molecular component is always small because of 
the small pre-factor l^ol/dj.. Even smaller P mo i is ex- 
pected along the state | b-f-c) because it has a weaker 
dependence on the x-axis as shown in Fig. [TJ Numeri- 
cally we find that the molecular population remains less 
than 1% for \r \/d± = 0.04. 

In a spherical harmonic trap, the entanglement prop- 
erties for ro — have already been studied before [lj|. 
Making use of our model as outlined above based on Ref. 
Q , we now extend the earlier result [l5[ to a cylindrical 
harmonic trap. To begin with, we first approximate the 
two fermion wave function ^ by neglecting the small 
molecular component, thus we obtain 



I*) = XN m < na mT a nll VaC > 



with the normalization constraint ^ r 



l'/n 



(10) 



= 1. 



r 7m,n ai 'e assumed real and symmetric for two identical 
atoms with rj m ll — r) n m . In the single atom basis state, 
the above wave function becomes 



I*) 



E 



77 m! „|min 2 ; 



I T1 I2) -II1T2) 
V2 



(11) 



Such a state has both spatial and spin degrees of free- 
dom, but the two degrees of freedom remain factorized. 
We therefore adopt the von Neumann entropy as our en- 
tanglement measure for a pure state. The factorized spin 
degree part simply contributes a ln(2), and the total en- 
tropy becomes E = E Bpat i a i+ln(2) [23|. While nontrivial, 
the spatial part of entanglement E spat i a i (abbreviated as 
pair entanglement) can be computed from the Schmidt 
decomposition (29j of the spatial wave function, i.e., we 
need to find 

^77m,n|m)i|n) 2 = K q|0q)l|^q>2, (12) 



For the two adiabatic states |a-d-b) and |b-f-c), we 
have numerically evaluated their spatial pair entangle- 
ment. The results are shown in Figs. [5] and [3] For both 
states, we find that pair entanglement is a smooth func- 
tion of —d±/a s . The pair entanglement first increases 
to some maximal value, then decreases and saturates to 
certain finite value. The overall dependence on the atom 
interaction strength remains essentially the same as be- 
fore [H. 
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FIG. 2: (Color online) Pair entanglement versus the atomic 
interaction strength at different trap aspect ratio A for the 
state |a-d-b) at |ro|/d± = 0.04. The dependence on ro is very 
small within the broad resonance regime. 
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with the Schmidt mode functions |0 q ) = \ipq), because 
Vm,n is symmetric. From this, the von Neumann entropy 
is found to be E spatial = - J2q K q m ( K q )- Computing the 
direct Schmidt decomposition for the three dimensional 
wave function turns out to be a quite demanding nu- 
merical task. Fortunately we can simplify this problem 
effectively to several one dimensional Schmidt decompo- 
sitions, which will be discussed elsewhere 30]. 



FIG. 3: (Color online) The same as in Fig. [2] except for state 
|b-f-c). 

We now discuss several interesting limits. First, 
for an approximately spherical trap with A ~ 1, the 
adiabatic state |a-d-b) corresponds to \a) or \b) at 
—d±/a s — > — oo or +oo, respectively. For A = 5/6, 
we find \b) « [E„=i, 2 (4J 2 ~ 2c t l2 4]|000) 1 |000) 2 , 
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whose pair entanglement is ln(2\/2) 
we again follow the notation of Ref. 



1.04. Here 

1 with 4. 

the creation operator for a fermionic atom indexed 
by n in the j-th trap direction. State |m.iTO 2 m3)„ 
therefore refers to atom n in the motional state 

= 7/6, we find 
] |000)i|000> 2 
1.39. For 
(c f ) 2 - 

-l,2;j—x,y 7 z \ nj J 



along the j-th. direction. For A 



|6>«[E„ = i, 2;J= ^(4-) 2 -2E 

with a pair entanglement In (4) 
A = 1, we find \b) oc [E, 



j=x,y c lj c 2j 



C lj C 2j 



1 000) 1 1 000) 2 with a pair entanglement 



ln(2-s/6) « 1.59. We find that the pair entanglement 
is always larger at the limit of a spherical trap with 
A = 1. For the adiabatic state |b-f-c), we find gen- 
erally that the pair entanglement at A = 1 is well 
separated from A ^ 1 because of the increased mo- 
tional state degeneracy. State |c) corresponds to the 
limit of —d±/a s — > +oo, which for A = 7/6 becomes 

|c) cc [E„=i, 2 (cD 2 - 2ct,4]|000)i|000> a , which is 
precisely the state \b) for A = 5/6 in the adiabatic 
state |a-d-b), whose pair entanglement therefore re- 



mains the same ln^v^) 
persists also for A = 5/6, 
same as state \b) for A = 



1.04. This correspondence 
where the state |c) is the 
7/6 in the adiabatic state 



|a-d-b). The state \c) at A = 1 is more complicatedly 
expressed as a linear combination of different single 
particle states as \c) oc [ - En=i.2:7=aM/.*( c L-) 4 



.t 



;nj) ( C nfe) 
t (J \3 , 

E« fc 4 7 -4 fc 4,-4]|ooo) 1 |ooo) 2 , 



-1.2; j—x,y.z \^nj ) 

i-) 2 (4^ 2 



®12j( C lj) \^2j) 
4 E« ;j yfc C lj C 2j ( C nfc) 

^ 1 i; ! i 2 / ' 2/,' wv,w / 1 w,n '/2* whose pair entan- 
glement is 551n(2) /24 + 71n(3) /8 - ln(5) /24 w 2.48. 

We next consider two extreme cases of A -C 1 or 
A ^> 1 corresponding to the quasi-one- and quasi-two- 
dimensional limits, respectively [22| . The pair entangle- 
ment results are shown in Fig. 2] for both A = 1/20 and 
A = 20. We find that the pair entanglement saturates 
to a lower value in the quasi-one-dimensional limit at 
A = 1/20, again because of the reduced motional state 
degeneracy. Our results here are of course limited to the 
validity regime of the model we adopt in the Hamilto- 
nian {!]). For a more rigorous treatment, please refer to 
the methods developed in Ref. [3l[ for harmonic traps or 
atomic waveguides. 



IV. A TOY MODEL 




FIG. 4: (Color online) Entropy as function of inverse scat- 
tering length for processes |a-d-b) and |b-f-c). Red solid lines 
are for A = 1/20, and blue dashed lines are for A = 20. The 
corresponding lower one is for process |a-d-b), while the upper 
one is for |b-f-c). 



For simplicity, we have truncated the motional states to 
only include the ground and the first excited state, de- 
noted respectively by |0) and |1) and separated by the 
trap frequency uj. The collisional interaction is mainly 
limited to the ground state manifold, giving rise to a level 
shift S when both atoms are in the ground state ( 1 00) ) 
and a simple atom excitation with strength r\. Double 
excitation is assumed small and neglected. 

This toy model can be easily solved, and the pair en- 
tanglement as function of interaction parameter tj/(lu — 5) 
is shown in Fig. [5] We see that the ground state 
is a simple product state of each atom in the mo- 
tional state |0) for -q = 0, i.e., displaying no entangle- 
ment. With the increase of atom-atom interaction, we 
find that the entanglement increases and saturates to 
21n(2) - ^ ln(2 + y/3) ~ 0.2458, where the approxima- 
tions for the toy become questionable. 

It is important to note that this toy model reproduces 
the same dependence of pair entanglement on the inter- 
action strength as for the system of two fermionic atoms 
during the BCS to BEC crossover. We thus feel it is im- 
portant to point out that from the point of view of two 
atom motional state entanglement, nothing particularly 
significant occurs during the BCS to BEC crossover. 



Much of our results above can in fact easily be appre- 
ciated from a toy model for two distinguishable atoms, 
as described by the Hamiltonian H = Hq + V, with Hq 
and V given by 

H = Rw(|10)ia(10| + |01}i2<01|), 
V = M|00)i 2 (00| 

+M|10) 12 (00| + |01)i2<00| + h.c). (13) 



V. CONCLUSION 

Before concluding, let us briefly summarize the gen- 
eralization of our current study to a narrow Feshbach 
resonance, where the bound state molecular component 
in Eq. @ cannot be neglected anymore. The pair entan- 
glement inside the molecular component can be included 
by performing an analogous symmetric Schmidt decom- 
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FIG. 5: (Color online) The dependence of pair entanglement 
as a function of interaction parameter t)/(lj — S). 

position on Eq. @. The total pair correlation can then 
be computed analogously in terms of the entropy of the 
independent Schmidt orbital expansion. Of course, such 
an approach would require more details about the model 
formulation and the exact determination of the molecular 
bound state wave function. 

In summary, we have studied pair entanglement be- 
tween two spin-1/2 fermionic atoms inside a single op- 



tical lattice site approximated by a cylindrical harmonic 
trap. We investigated thoroughly the dependence of pair 
entanglement on the trap strength and geometry and on 
atom-atom interaction strength along the complete BCS 
to BEC crossover and focused on the broad resonance 
regime. We developed a formalism for studying pair en- 
tanglement including the effect of an effective range rg 
for two interacting atoms at low energy. In the limit of 
a broad Feshbach resonance, where the effect of tq be- 
comes negligibly small, our result reduces to the theory 
developed before for evaluating pair entanglement in a 
single open channel of two atoms without the presence of 
a bound molecule state. We find that pair entanglement 
changes significantly against atomic pair interaction be- 
cause of the induced motional orbital deformations. In 
general, however, the exact value of the spatial pair en- 
tanglement is also governed by the motional state de- 
generacies. As a rule of thumb, we find that spherical 
harmonic traps generally give rise to larger pair entan- 
glement. We hope our study will provide new insights 
into the applications of quantum degenerate lattice sys- 
tems to quantum information science. 
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